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The IceCube detector allows for the first time a measurement of atmospheric muon and neutrino energy spectra 
from tens of GeV up to the PeV range. The lepton flux in the highest energy region depends on both the primary 
cosmic ray composition around the "knee" and the contribution from prompt decays of mostly charmed hadrons 
produced in air showers. It is demonstrated here that a direct measurement of the atmospheric muon spectrum 
in the region above 100 TeV is feasible using data that is already available. 



1. INTRODUCTION 

Neutrino astronomy is a relatively young 
branch of astrophysics. During the past two 
decades, several prototype detectors were built 
(Baikal, AMANDA, ANTARES). This estab- 
lished the basic techniques of neutrino detection 
in the Very High Energy (VHE) range from tens 
of GeV to PeV and above, but did not lead to 
any discovery of an astrophysical signal. 

IceCube is the first of a new generation of neu- 
trino telescopes that will have instrumented vol- 
umes of the order of one cubic kilometer. Once 
completed, it is expected to put detection of as- 
trophysical point sorces within reach [1 . There 
are three distinct components in the IceCube ar- 
ray. The main volume detector Inlce consists of 
80 strings each holding 60 Digital Optical Mod- 
ules (DOMs), located in the ice at a depth of 1.5 
to 2.5 km. The surface air shower array IceTop 
is comprised of 160 surface tanks, and six more 
strings with closer-spaced DOMs will be deployed 
as the low-neutrino-energy extension Deep Core. 
A schematic drawing of the detector is shown in 
figure [l] As of mid- January 2009, almost sixty 
strings and a similar number of surface stations 
have been deployed. The entire system is ex- 
pected to be complete after the austral summer 
2010/11 at the latest. 

The main channel for neutrino detection is the 
process ^ N /i + X. The muon in the fi- 
nal state emits Cherenkov radiation when passing 




Figure 1. Schematic view of the IceCube detector. 
The Inlce neutrino detector is located from 1.5 to 
2.5 km under the surface of the ice. For top and 
bottom of the detector, angles corresponding to 
slant depths of 10 km and 20 km are indicated. 
The region in between holds the greatest potential 
for a direct atmospheric muon energy spectrum 
measurement. 



* www. icecube, wise. edu/collaboration/authorlists/2008/4. html 
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Figure 2. True angular distribution of muon 
tracks from cosmic-ray induced air showers and 
atmospheric neutrinos for trigger-level Monte 
Carlo. 



through ice, which can be measured using photo- 
multiplier tubes. The direction of the muon track 
can then be reconstructed by means of a fit to the 
arrival time of the Cherenkov photons. Also, in- 
tegration over the measured charge yields calori- 
metric information that can be used to estimate 
the muon energy and by extension that of the 
parent neutrino. 

By far the dominant background for neutrino- 
induced muons comes from cosmic-ray air show- 
ers. Figure [2] shows the respective event rates 
from atmospheric muons and atmospheric neu- 
trinos. Above the horizon, muons dominate by 
several orders of magnitude. Extraction of a neu- 
trino signal, other than at the very highest ener- 
gies (tens of PeV and above), is therefore effec- 
tively only possible below the horizon, where the 
Earth itself can be used as a filter to attenuate 
cosmic-ray induced muons. 

Even below the horizon, the background from 
misreconstructed downgoing muons is several or- 
ders of magnitude higher than the neutrino event 
rate. Figure |3] illustrates the experimental sit- 
uation. The main challenge for any VHE neu- 
trino detector is to find a set of quality cuts that 
reduces the amount of misreconstructed muon 
tracks by four or five orders of magnitude while 



Figure 3. Reconstructed muon track direction in 
IceCube at trigger level, before any quality cuts. 



preserving as much as possible of the neutrino 
signal. 

2. MUON DATA IN ICECUBE 

However, the large number of atmospheric 
muon events registered by the detector should not 
simply be seen as an undesirable background for 
neutrino searches. In fact, it represents a high- 
statistics data set that can be used for calibration 
and verification of the detector system as well as 
physics investigations. 

IceCube data are transmitted to the northern 
hemisphere in near-real time in order to allow im- 
mediate high-level processing. Since all commu- 
nication with the South Pole takes place over Low 
Earth Orbit satellites, the maximum data rate is 
restricted to about 32.5 GB/day, corresponding 
to roughly 6% of the total event rate of 1400 Hz in 
the 40-string configuration. Data reduction is ac- 
complished using filter conditions that make use 
of simple quality parameters from online event re- 
constructions. Each condition leads to a distinct 
filter stream. For atmospheric muons, the rele- 
vant data sets are: 

• Muon Filter: This stream contains all events 
reconstructed near or below the horizon. Its 
main purpose is the analysis of up-going 
neutrino-induced muon tracks. Event rate: 
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IceCube-22 Data vs. Monte Carlo Simulation Data | 
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Figure 4. All- Sky Muon Flux as measured with 
IC-22, renormalized to trigger level (actual event 
rates after cuts are 8%-20% of the value shown). 
The data in this plot represent a total livetime 
of 275.6 days. The downgoing muon flux has 
been simulated using CORSIKA/SIBYLL and is 
shown separately for single and multiple-primary 
(coincident) cosmic ray showers. Muon neutrino 
events were simulated up to cos Ozen = 0.25. 



tens of Hz 

• Minimum Bias: For general quality con- 
trol and calibration, an unbiased sample of 
events is provided. In the 22-string configu- 
ration (IC-22), every 100th event was stored. 
This has since been reduced to every 2000th. 
Event rate: < 1 Hz 

• nanoDST: For every trigger, a rudimentary 
set of low- level event parameters is kept. Ex- 
amples are the total number of photoelectrons 
registered and track parameters from a basic- 
level reconstruction. Event rate: > 1 kHz 

• IceTop Coincident: This stream contains 
Inlce events that are coincident with air show- 
ers triggering the IceTop surface array. Due 
to the placement of IceTop, this data sam- 
ple is restricted to the area around the zenith. 
Event rate: tens of Hz 

Using a combination of minimum bias and 
muon-filtered data, the performance of IceCube 
has been verified by comparing the measured 
muon track zenith angle distribution to a Monte 



Carlo simulation. Figure [4] shows the result, rep- 
resenting the first time in IceCube/ AMANDA 
that the muon fiux was measured over the en- 
tire sky using a uniform set of quality parame- 
ters. Data and simulation generally agree to bet- 
ter than 20%, meaning that detector systematics 
are known reasonably well. A similar measure- 
ment had earlier been performed for AMANDA 
[2j, but was restricted to angles above the hori- 
zon. 

The high statistics available in the nanoDST 
format allows to search for possible anisotropics 
in the primary fiux. Both the Tibet Air Shower 
Array and MILAGRO report a dipole anisotropy 
at primary energies in the TeV range |3|4j . Ice- 
Cube should be able to confirm this effect based 
on data taken in its 22-string configuration, which 
took data from May 2006 to March 2007. An in- 
vestigation into this issue is ongoing. 

The main goal of IceTop is to measure the com- 
position of cosmic rays around the "knee", lo- 
cated near Eprim = 3PeV. An important param- 
eter here is the ratio of electromagnetic to muonic 
particles in an air shower. This ratio can be esti- 
mated by combining information from the surface 
array and the Inlce detector. A preliminary mea- 
surement of the primary cosmic ray spectrum has 
recently been completed, confirming that IceTop 
performs according to expectations 

3. MUON SPECTRUM MEASURE- 
MENT 

3.1. Motivation 

The origin of the "knee" in the cosmic spec- 
trum can be explained if the spectrum of the in- 
cident cosmic ray primaries depends on the en- 
ergy per nucleon, with a sharp steepening at 
Ecr/A 4.5 PeV [6 . This change in compos- 
tion is consistent with measurements by the KAS- 
C ADE array [7 . Since the ratio of median parent 
cosmic ray and muon energy is < 10 at VHE en- 
ergies [8^ , a steepening of the energy spectrum by 
nucleon of cosmic rays at a few PeV will have a 
measurable effect on the atmospheric muon spec- 
trum at energies of hundreds of TeV. 

It has also been pointed out that the muon fiux 
measured with neutrino telescopes will be sen- 
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sitive to the prompt flux from charm decay in 
air showers [9 . The ratio of prompt neutrinos 
to prompt muons is close to unity, whereas con- 
ventional muons outnumber neutrinos by a factor 
of 3-5 at TeV energies, making detection of the 
prompt component more difficult for muons than 
for neutrinos. A review of prompt lepton flux cal- 
culations is given in [10 . AMANDA already set a 
limit on prompt neutrinos [11 . Since the prompt 
contribution to the atmospheric flux represents 
an important background in searches for a diffuse 
astrophysical neutrino flux, an independent limit 
from muons would be of particular importance to 
neutrino telescopes. 

In order to measure any deviations from ex- 
pectation, it is essential to accurately predict 
the spectrum of muons from light meson de- 
cays. While muon neutrinos at TeV energies 
mostly come from the process K + X, 

for kinematical reasons muons originate predom- 
inantly in pion decays n ^ ^ ji [12]. An es- 
timate of the pion production cross section from 
accelerator experiments gives an uncertainty of 
15% + 12.2% • /o^io(^,r/500 GeV) at xiah > 0.1 
above 500 GeV [T3^. This value should also ap- 
ply in good approximation to the conventional 
(non-prompt) muon flux. The relatively small 
uncertainty for the conventional flux can be used 
to calibrate any detector response to high-energy 
muons. 

The energy spectrum of atmospheric muons has 
never been measured beyond a few tens of TeV. 
All results so far are consistent with, or slightly 
above, theoretical expectation [14 . As will be 
shown, IceCube should be able to extend the en- 
ergy range for muon spectrum measurements by 
at least an order of magnitude. 

3.2. Methodology 

To be registered by Inlce, muons have to tra- 
verse at least 1.5 kilometers of ice and still be 
energetic enough to trigger the detector. This 
means that there is an effective energy threshold 
in IceCube requiring a minimum muon energy of 
about 400 GeV for vertical air showers. At larger 
zenith angles and increasing depth, this threshold 
energy increases correspondingly [TS]. 

Since the granularity of high-energy neu- 
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Figure 5. Muon multiplicity averaged for all cos- 
mic ray air showers at closest approach to the 
center of the Inlce detector for three different 
slant depths. All figures were derived from a 
CORSIK A/SIB YLL-based Monte-Carlo simula- 
tion and apply to IC-22. The primary cosmic 
ray spectrum was parameterized according to the 
poly-gonato model [6]. Propagation of muons 
through ice was simulated using MMC [15]. 



trino telescopes prevents resolution of individual 
muons, only the total energy of a shower can be 
measured. This fact has been used to derive the 
spectrum of the cosmic ray flux using AMANDA 
muon data [16]. However, as shown in figure 
[5) the muon multiplicity decreases substantially 
with increasing slant depth. At a depth of 10 km 
water equivalent (kmwe), 90% of all triggering air 
showers only contain a single muon at closest ap- 
proach to the center of the Inlce detector. 

After a sufficiently strong quality cut, the de- 
pendence of the average muon surface energy on 
the measured slant depth agrees well with the 
true value, as shown in fig. [6j The deep edge 
of the region useful for a muon spectrum mea- 
surement is determined by two limiting factors. 
First, since the slant depth at a given zenith an- 
gle 6zen is proportional to l/cos6zen^ finite angu- 
lar resolution {SO ~ 0.7° in IC-22) in combination 
with steeply decreasing flux rates will always re- 
sult in a point where the muon event rate becomes 
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Figure 6. Average surface energy of the most 
energetic muon in air showers depending on the 
slant depth measured from the surface to the 
point of closest approach to the center of the Inlce 
detector. 
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Figure 8. Total IC-22 event yield from cosmic-ray 
induced atmospheric showers above given surface 
energy of the most energetic muon. The three 
different curves represent all IceCube triggers, 
events included in the muon filter event stream 
and after application of stringent quality cuts. 
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Figure 7. Energy of most energetic muon in air 
showers at the point of closest approach to the 
center of the Inlce detector. The energy resolu- 
tion of Inlce has been estimated as 0.3 (0.5) in 
logioE for = 10 TeV (1 TeV) 



dominated by tracks from higher angles. This is 
illustrated by the turnover in the curve represent- 
ing the measured depth after quality cuts. Sec- 
ond, atmospheric neutrinos become an increas- 
ingly important background, eventually drowning 
out the muon signal. Taken together, these two 
factors lead to an upper bound for the useful slant 
depth region of slightly less than 20 kmwe. 

For an accurate energy measurement it is also 
necessary for the muons to retain a significant 
fraction of their surface energy. Muon ener- 
gies can effectively only be determined for ener- 
gies above 1 TeV, where the energy loss dE^/dx 
strongly depends on E^ Figure [t] shows that 
even after passing through >10 km of ice, muons 
with Esurf > 100 TeV still fulfin this condition. 
Given sufficient statistics, a measurement of slant 
depth and energy therefore provides sufficient in- 
formation to derive the muon energy spectrum. 

Figure |8] shows the integrated event yield above 
a given muon surface energy for three different cut 
levels in IC-22. Even at this early stage in the 
construction of the detector, sufficient statistics 
should be available up to muon energies of hun- 
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dreds of TeV. Since the measurement takes place 
at angles near the horizon, the conventional muon 
flux is maximized. This favors an investigation 
into the cosmic ray composition over a search for 
prompt muons. 

It was further proposed to use topological sep- 
aration of muon tracks belonging to the same air 
shower to investigate high.-pt{> 3 GeV/c) muon 
tracks in IceCube [18 . This data can be used for 
an independent measurement of the cosmic ray 
composition which should mutually validate the 
result from the analysis oulined above. 

4. SUMMARY AND OUTLOOK 

Far from simply constituting a nuisance in the 
detection of neutrinos, atmospheric muons pro- 
vide a unique opportunity to extend the scientific 
value of large- volume astrophysical neutrino tele- 
scopes. The combination of km^-scale detector 
volume and large amount of shielding material 
provides an opportunity to directly measure the 
atmospheric muon spectrum at low zenith angles 
in the range from tens of TeV up to 1 PeV and 
beyond. 

Detection of a steepening of the muon spectrum 
in this energy region would provide independent 
evidence for a change in cosmic ray composition 
that does not rely on knowledge of the height of 
the shower maximum in the atmosphere, as in 
air shower arrays. Also, a limit on prompt muon 
fiuxes, and hence on charm production in nuclear 
collisions at Eqm ^ 1 TeV can be set. At these 
energies, accelerator measurements are only pos- 
sible in collisions, and thus restricted to high pt 
values [19 . Neutrino detectors allow investiga- 
tion of the forward region and provide a natural 
complement to laboratory experiments. 

Upcoming fixed-target experiments at the LHC 
will investigate light meson production at TeV 
energies [20]. In combination with results 
from second-generation neutrino telescopes, this 
should have a considerable impact on a vari- 
ety of astrophysical investigations. These range 
from diffuse astrophysical neutrino searches to 
the measurement of cosmic ray fiuxes at UHE 
energies with large air shower arrays. Further 
opportunities for physics investigations will un- 



doubtedly arise as neutrino detectors increase in 
size and the field becomes more mature. 

The analysis of IceCube data with the goal of 
deriving an atmospheric muon spectrum is ongo- 
ing and should yield results in the near future. 
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